Ferromagnetic transition in a double-exchange model 
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We calculate the temperature of a ferromagnetic transition in a double-exchange model with 
classical core spins for arbitrary relation between Hund exchange coupling and electron band width 
by solving the Dynamical Mean Field Approximation equations. 
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I. INTRODUCTION 

The double-exchange (DE) model is one of the 

basic ones in the theory of magnetism. Magnetic ordering 
appears in this model due to Hund exchange coupling be- 
tween the core spins and the mobile carriers. The Hamil- 
tonian of the model is 
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where c and are the electrons annihilation and cre- 
ation operators, S„ is the operator of core spin, is 
the electron hopping, J is Hund exchange coupling be- 
tween a core spin and a conduction electron, a is the 
vector of the Pauli matrices, and a, (3 are spin indices. 

The model (the core spins being treated as classical 
vectors) was thoroughly studied already in the papers 
During the last years, because of a general inter- 
est in manganites, the model was brought in the focus of 
attention, and a lot more was achieved (see reviews j§,|]]6) 
and references therein). However, some basic properties 
of the model are still known only partially. For example, 
most of the papers dealing with the DE model, starting 
from classical paper by De Gennes ||], considered the 
DE Hamiltonian with infinite exchange (and with the 
addition of the antiferromagnetic superexchange, which 
is crucial for the explanation of magnetic properties of 
manganites). 

The other extremity which was studied is the partic- 
ular case of weak exchange (much less than the electron 
bandwidth), when the DE Hamiltonian can be reduced to 
Ruderman-Kittel-Kasuya-Yosida (RKKY) Hamiltonian 
(see review Q and references therein). 

In this paper we calculate the temperature of a ferro- 
magnet - paramagnet transition T c in a double-exchange 
model for arbitrary relation between Hund exchange cou- 
pling and electron band width by solving the Dynami- 
cal Mean Field Approximation equations. Note that we 
treat the core spins as classical vectors. (When the quan- 
tum nature of the core spins is taken into account, the 



Hamiltonian ([[]) , which in this case is often called the pe- 
riodic Kondo model, becomes much more complicated; 
only scanty results were obtained for the model up to 
now.) The problem of classical spins, as we shall see, 
combines tractability with rich and interesting physics. 



II. HAMILTONIAN AND DMFA EQUATIONS 

Like it was said above, we consider spins as classical 
vectors S„ = m„ with the normalization |m| 2 = 1, Thus 
the DE Hamiltonian in a single electron representation 
can be presented as 
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We have a problem of electron scattered by core spins, 
the probability of any given core spin configuration de- 
pending upon the energy of electron subsystem. To solve 
the problem we will use the Dynamical Mean Field Ap- 
proximation (DMFA) (see |||| an d references therein). 

In this approach, first we calculate an averaged, with 
respect to random orientation of core spins, density of 
states of electron in a random core spins configuration, 
treating electron scattering in a single site approxima- 
tion, and considering the probability of any configuration 
as given. We introduce Green's function 



G(E) = (E-H)-\ 
In this approximation the averaged locator 

G\ oc (E) — (G nn (E) 
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(4) 



is expressed through the the local self-energy S by the 
equation 



where 



G loc (£)=So (e-E{E) 
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is the bare (in the absence of the exchange interaction) 
locator. The self-energy satisfies equation 
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where (X(m)) = J X(m)P(m), and -P(m) is a probabil- 
ity of a given spin orientation (one-site probability). The 
quantities G and £ are 2x2 matrices in spin space. 

In PM phase P(m) = const, the averaging in Eq. 
(0) can be performed explicitly, S = Ei, G = gl = 
go(E — 11)1, where 1 is a unity matrix, and we obtain 
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To first approximation of the DMFA, leading to Eq. 
(0) , has a simple physical meaning. We reduce the prob- 
lem of electron scattering due to many spins, each with 
the scattering potential Jm ■ a, to a problem of a scat- 
tering due to a single spin with the effective scattering 
potential — Jm -a— S, embedded in an effective medium, 
described by the Hamiltonian tk + S, and, hence, by lo- 
cator G\ oc . 

The same MF approach leads to the second DMFA 
approximation - the approximation for the one-site prob- 
ability -P(m), which allows to perform averaging in FM 
phase. Consider again a single spin with the effective 
scattering potential in an effective medium. The change 
in the number of states of the electron gas due to such 
spin is 
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where Y + = Y(E + iO). So the change in thermodynamic 
potential is 0HH 



(10) 



/3Afi(m) = J f(E)AN(E,m)dE, 



where f(E) is the Fermi function, the chemical potential 
is found from the equation 



2 f°° 
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and n is the number of electrons per site. The result for 
the one-site probability reads: 



P(m) tx exp [-j3Afi(m)] 



(12) 



Eqs. ([?]) and (|T^) are the system of two non-linear (in- 
tegral) equations for S(£') and P(m), which one should 
solve to find thermodynamic properties of the model. 
However, in linear approximation with respect to macro- 
scopic magnetization M, we can reduce this complicated 
system to a traditional MF equation for M p5[: 



P(m) cx exp (-3/3T c M • m) 



(13) 



The parameter T c is formally introduced as a coefficient 
in the linear term of the expansion of Af2(m) with re- 
spect to M (the reason for the notation we have chosen 
and for the numerical coefficient 3 will be clear immedi- 
ately); it is determined by the properties of the system 
in paramagnetic phase. Non-trivial solution of the MF 
equation 



M = (m) 



(14) 



can exist only for T < T c , hence T c the Curie tempera- 
ture. 



III. T c FOR SEMI-CIRCULAR DOS 

For simplicity consider the semi-circular (SC) bare den- 
sity of states (DOS) N a (e), the bandwidth being 2W. 
Then 



9o(E) = 



N (e)de 



E-e 
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(15) 



For this case 



£(E) = E - 2wG loc - G- a l, (16) 

where w = W 2 /8, and Eq. (0) and (g) take respectively 
the form 
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Expanding Eq. ( |l7| ) and then Eq. (||) with respect to 
M, after straightforward algebra, for the T c we obtain 



T r . = 
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Formally speaking because the integral in Eq. (19) con- 
tains Fermi function critical temperature enters also into 
the r.h.s. of the equation. But in all cases T c turns out to 
be much less the chemical potential, so we can consider 
electron gas as degenerate, and Eq. ([l9]) is an explicit 
formula for calculation of T c . 

Eq. (^|) is the main result of our paper. The anal- 
ysis of this equation let us start from the limiting case 
J = oo. In this case integral can be calculated explicitly 
and we obtain |T4l 
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where y is an implicit function of concentration, given by 

1 _ l 

2 vr 
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The result 



equation n = | — ± (^sin 1 y + yyl- 

coincides with that known previously |9[| . 

For arbitrary exchange the integral can be calculated 
only numerically, but before we present the results of 
calculations we should state that in part of the J/W — n 
plane Eq. (|l|) gives T c < 0. In fact, like in any MF 
theory of a second order phase transition, in our calcu- 
lation of the critical temperature we started from a high 
temperature paramagnetic phase and decreasing the tem- 
perature were looking for an instability of a model with 
respect to appearance of small spontaneous magnetic mo- 
ment that is for the appearance of a nontrivial solution 
of the MF equation. Negative T c in some part of the 
J/W — n plane means, that at any temperature, includ- 
ing T — 0, the paramagnetic phase here is stable with 
respect to the appearance of small spontaneous magnetic 
moment even at T = 0, and excludes ferromagnetism in 
that region of the plane. The part of the critical sur- 
face corresponding to the region of the parameters plane 
where T c > is presented on FIG. 1. 
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FIG. 1. T c as a function of relative strength of the Hund 
exchange J/W and electron concentration n. 



In fact, this equation is the MF approximation plj] for 
the RKKY Hamiltonian M , to which the original Hamil- 
tonia n (|l| ) can be reduced to in the case J <C W. So 
Eq. (|2l|) docs not involve either the coherent potential 
approximation (Eq. (0)), or the SC density of states. 

Anyhow, for the SC density of states we use, Eq. (|2l| ) 
gives ferromagnetic ground state for n < 0.4, which qual- 
itatively agrees with the result of numerical calculations, 
giving FM ground state for n < 0.25 for the three prin- 
cipal cubic lattices fll6|| . 

To formulate the second argument, let us compare the 
energies of the PM state 



E„ 



E Im g + dE; 



(22) 



and of the saturated FM state 



E f = I ' F E ^2n (E±J) dE, 
J-oo (±) 



(23) 



where Ep is the appropriate Fermi energy. On fig. 2 we 
plotted simultaneously the curve given by by equation 
E p — Ef and by equation T c — 0, where T c is obtained 
fromEq. (ph. 
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FIG. 2. The FM region boundary in the coordinates of rel- 
ative strength of the Hund exchange J/W and electron con- 
centration n. 



IV. DISCUSSION 

Let us finally discuss, whether the PM-FM transition 
observed with decreasing temperature can be preceded 
by the transition from the the PM phase to some mag- 
netic phase other than FM (say, antiferromagnetic)? We 
would like to present some heuristic arguments that this 
is not the case. 

First, consider the case of weak exchange J W . In 
this case Eq. (n9h takes the form 



2 J 2 f°° 
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— — Im g (# 4 
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dE. (21) 



The close vicinity of the abovementioned curves sup- 
ports the belief that the curve T c = is the quantum crit- 
ical line (see JItJ and references therein), which bounds 
ferromagnetic phase. Also, the boundary of ferromag- 
netic phase on FIG. 2 agrees with those obtained on the 
basis of numerical calculations and from qualitative 
reasoning |p^| . 

The destruction of ferromagnetic ground state occurs 
because finite double exchange between the itinerant elec- 
trons and core spins, unlike an infinite one, by itself gen- 
erates effective antiferromagnetic exchange between the 
core spins (which was absent in original Hamiltonian). 
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Due to the fact that our main result (equation for the 
transition temperature) indicates it's own limits of valid- 
ity, we can find the boundaries of ferromagnetic phase 
without analyzing what phases are beyond the bound- 
aries. 

Finally, we would like to mention, that in the DMFA, 
as one can easily see from Eq. (^), density of electron 
states in paramagnetic phase does not depend upon elec- 
tron concentration. In this case, the derivative of chem- 
ical potential with respect to number of electrons is just 
the inverse density of states at the Fermi level (for the 
degenerate electron gas), and is always positive. Hence, 
there is no phase separation in paramagnetic phase. 

In conclusion, we explicitly formulated the Dynami- 
cal Mean Field Approximation equations for the dou- 
ble exchange model with classical spins for arbitrary re- 
lation between Hund exchange and the electron band- 
width. Near paramagnetic-ferromagnetic transition crit- 
ical point, these equations were reduced to a MF equa- 
tion, describing a single spin in an effective field, propor- 
tional to the macroscopic magnetization. The effective 
exchange interaction, entering into the MF equation was 
found for the semicircular electron density of states. We 
thus calculated the transition temperature T c as a func- 
tion of Hund exchange interaction and electron density 
in the whole parameters plane. The results obtained also 
allow to plot the boundaries of the ferromagnetic region 
on the model phase diagram. 
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